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ON THE COVER — The area of the heavens around the Orion Constellation, 
shown in the cover photograph made through the 120-inch telescope of the Lick 
observatory, is also the region of observations with an infrared telescope de- 
veloped by University of Minnesota astro-physicists. The infrared sensory equip- 
ment reveals stellar bodies that could not be studied by conventional telescopes, 
and it is expected to provide data on the birth of stars. 
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FOREWORD 


This report covers the first twelve-month period, 

1 July 1972 to 30 June 1973 of support under NASA Grant 
NGL 24-005-263. The purpose of this project is to foster 
and develop new applications of remote sensing under an 
interdisciplinary effort with the general title "A Study of 
Minnesota Forests and Lakes Using Data from Earth Resources 
Technology Satellites . 11 The following five reports reflect 
the progress recently reported prior to this date and the 
work proposed by new investigators for the coming year. 
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" Remote Sensing in Lake Superior Studies . ” 

Investigator: Dr. Michael Sydor , 

Department of Physics ( 

The surface water transport in the extreme western ann of Lake Superior 

is characterized by a localized surface current loop, and a general eastward 

♦ 

transport of the surface water out of the arm along the south shore of the 
lake.^ This is accompanied by a cold water upvelling along the North Shore 
of Lake Superior. Our activities in the summer of 1971 and 1972 supported 
under a NASA Remote Sensing Program were concerned with the investigation of 
the transport processes at the tip of the lake in greater detail. Our project 
is a part of the continuing Lake Superior Studies at the University of Minne- 
sota. The results of our water current measurements, and the measurements on 
thermocline activity at three stations located along the northern shore of 
Lake Superior from Knife River to Duluth, show that the bottom currents have 
a prevailing northerly direction at station 1, Fig. 1, displaying a broad area 
of upvelling south of Knife River. A veil defined current pattern, somewhat 
smaller than anticipated from previous studies by T. Olson and T. Odlaug'*' 
develop es in the area between the upvelling and Duluth. The upvelling serves 
as a turn around path for the cell currents. The current pattern was investi- 
gated using water current and temperature profile measurements in conjunction 
with ERTS Satellite data. It was found that an eddy loop developes near Duluth 
because of warn water pileup at the tip of the lake. The circulation pattern 
associated with the warm water cell displays a characteristic feature influenced 
by the winds and the depth of the cell. The circulation persists until the warm 
water is driven out of the area. 


1. Lake Superior Studies, T. Olson et al, Minnesota Public School of Health. 

2. Charles E. Adams, Jr., Proc. 13 Conf. Great Lakes Res., (862-879) 1970. 
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The knowledge of the summertime current patterns is important in that 

I 

these currents are responsible for the sweeping out of the polluted waters 
near* Duluth. The water currents tend to disappear in the winter, particularly 
with the development of ice cover on the lake, causing accumulation of polluted 
waters from the St. Louis Biver at the tip of the lake. This situation is 
potentially hazardous since three water intakes are located in the area. 

Little information, even of qualitative nature, was available on the winter 
water transport in the lake. The winter currents were very low and were dif- 
ficult to measure on account of the lack of safe ice cover. Thus other methods 
for determining the winter transport had to be established. Since the water 
quality parameters for the St* Louis River are considerably different from 
those of Lake Superior, and the water quality at the Cloquet water intake 
located 2 miles off shore at the tip of the lake, could be monitored on a 
continuous basis, we decided to use the properties of the river water to i- 

dentify suitable parameters for measurements of the rate of accumulation of 

\ 

the effluents under the ice. For this purpose, we- had to examine the prop- 
erties of the river water and its transport in the bay area first. Such data 
for the winter months was . unfortunately nonexistent, thus we had to extend 
our measurements to the Duluth-Superior harbor area. The results of our 
studies give the yearly variation in water quality parameters for St.. Louis 
River, indicating that the electric conductivity, transparency and dissolved 
oxygen could well serve as effluent tracers. These parameters were subse- 
quently used to measure the effluent accumulation at the Cloquet water in- 
take during the time of ice cover on the lake. Some results for measurement 
of ice cover are also discussed, including a prediction for ice breakup in 
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the harbor. We vrere also involved in evaluation of an- experimental bubbler 
system installed in December 1972 near the Superior entry to the lake. Only 
some of the data for the past year is presented here, sufficient ve hope to 
serve as a progress report justifying our support from NASA for this program, 
and our subsequent request for manned aircraft remote sensing data for our 
area. Consideration for the practical aspects of our investigation to the 
community is outlined as required by the contract. 
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Results end Discussion 

A series of temperature and vater velocity profile measurements were 
made from July 25 to September- 15, 1972, at the three stations located in 
the extreme western arm of Lake Superior to determine the nature of currents 
at the tip of the lake. Temperatures were measured using pressure operated 
bathythermograph usable down to 200 1 and thermistor probes for entire temp- 
erature profile. Water currents were measured using vand type current meters. 

From the end of July until August 9, no pronounced steady stratification 
took place. A series of days vith variable vinds produced, under the in- 
fluence of N.E. vinds, varm vater pileup at the tip of the lake, which vas 
subsequently driven out by westerly vinds. Considerable mixing took place 
in that period of time as shown by the temperature profiles taken at station 
01 off Knife River for August 1-7, Fig. 2. Cold vater near the shore vas 
often closer to the surface indicating presence of upvelling. The currents 
below iJhermocline at station #1 often pointed due north at speeds V2 cm/sec. 
These seemed to be characteristically present at Station 01 through August 
and September. The upvelling diminished for times of prolonged N.E. vinds. 
This broad area of upvelling off Knife Paver tends to contain the warm 
vater cell developed at the tip of the lake. Such vas the case for August 
9 data when a distinct stable stratification developed after the prolonged 
N.E. vinds for August 6-8. The temperature profile for the day Fig. 2, . 
shows the stratification vith slight upvelling detected by comparing the 
temperature profile for points 1 and 3 miles off shore. The accompanying 
thermocline activity Fig. 3 vas characterized by internal waves which had 
a rather regular ^ 50 min. period. Fig. 3. The crest of the wave vas calcu- 
lated to be ^ 10 ft. , using the temperature profile taken at the time of the 
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crest of internal wave and at the slack time. The water velocity at depth 
below the thermocline was ^ 1 cm/sec. due'N.W. The surface water velocity 
was towards vest at ^ ^ cm/sec. 

The BRTS image for August 12 shows the turbidity pattern for the period. 
The turbidity resulted from the lake bottom activity developed during the 
N.E. winds. The turbidity pattern indicates an eddy current at the very tip 
of the lake as shown in Fig. 4 and 5* This type of current pattern prevailed 
in the area from then on into August and September, and was verified with the 
current measurements at the stations. The stratification seen on August 9, 
persisted for August 11 though on that day the thermocline activity (Fig. 3) 
was not characterized by regular internal surges , but showed a steady shift 
in thermocline, as the cell deepened. Subsequent easterly winds caused an 
even further stratification as shown by temperature records from August 11 
to August 23. Considerable precipitation occurred in mid-August, causing 
extensive turbidity plumes from the rivers. The edge of the turbidity plume 
off Knife Hirer, on August 17>.was .7 of a mile away from the shore. The 

transparency profile of plume, measured using a homemade transparency meter 
(20 cm. light path, broad sensitivity peaking at 6000 X. } shows that the 
plume edge consists of a tongue of ^10 ft. layer of turbid water riding 
up on top of the upwelling water. Fig. 5. (Secchi disc transparency on the 
clear side of plume was 25 ft., while inside the plume it was ,2.5 ft. The 
plume deepened as the shore was approached, presumably from 'settling of 
coarser material. A presence of a separate shallow warm water layer is 
seen in temperature profile at Station Hi for August IT and for August 23, 
where it is still discernible but somewhat obscured by a much deeper layer 
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of warm vater which vas driven in "by prolonged easterly vinds following a 
torrential rain for August 20. The edge of this varm water bounded the 
region of the up welling and displayed long Langmuir cells running parallel 
to the shore for one to two miles. The winds were consistently easterly 
until August 27. The profile for August 28 (Fig. 6), shows the strati- 
fication at the end of the month. ERTS data for August 30 again shows a 
well defined eddy loop quite similar to the August 12 pattern. The winds 
were westerly for the day, and the up welling along the North Shore vas 
pronounced as shown by the temperature profiles. . Water velocity profiles 
are shown for August 30 and September 1. ERTS records of the turbidity 
pattern confirmed the direction of the measured surface flow. A comparison 
of temperature for a given station on these two days shows that the cell of 
warm water vas still confined to the tip of the lake by the upvelling. The 
eddy current persisted further into the season in spite of westerly winds, 
but the cell became shal lover as the warm water was transported out of the 
area as indicated by temperatures for September 7-12. The eddy current 
vas again evident in the October 6 ERTS image, due to the prolonged easterly 
vinds late in September and early October. The November 29 ERTS image shows 
that the loop current disappears altogether in agreement with the wind data 
which showed h days of strong westerly vinds prior to the day the image vas 
taken. No temperature profile records are available past September lU. 

s ■ 

The data to date, tentatively established the general nature of the 
summer-fall currents at the tip of the lake. More stations are needed to 
measure the vater velocity at sufficient number of points necessary for an 
analytical treatment of the transport. The proposed stations for future 
studies are indicated in Fig. 1. 
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To study the wintertime water transport in the area, through tracing 

• - * * f 

M f 

of effluents, it was first necessary to establish which parameters were 
suitable as tracers. It was also necessary to examine the transport process 
vithin the Superior Bay itself, as the water quality parameters in the bay 
are influenced by lake water surging into the bay through the entries. 

Chemical analysis of samples would be the most reliable method for tracing 
the effluents. Unfortunately it quickly defeats any attempt at extensive 
measurements since the chemical analysis is tedious. We thus spent our effort 
on identifying the parameters which could be suitable for relatively quick 
measurements of concentration of the effluents in the lake. Figures 12-19 
show the results for some of the water quality parameters measured on weekly 
basis for the year 1972. They indicate the expected correlation between 
suspended solids and turbidity, both showing the spring runoff effects. The 
temperature measurements indicate considerable inflow of lake water into the 
Duluth . harbor area in May and June, Figure 1 6. The N.E. winds kept the ice 
cover at the tip of the lake until June 9> 1972. - The turbidity generally 
dropped for the winter months when considerable settling occurred. Pro- 
nounced stratification with depth occurred as a result of settling of pol- 
lutants in areas where industrial dumping takes place. Drop in the dissolved 
oxygen at the bottom of the shipping channel was much greater than indicated 
in the seasonal variation shown in Fig. 15, which pertains to water values at 
5 ft. depth. Table 1 in Figure 20 shows an example of chemical analysis for 
the surface and the bottom water in the channel. Figure 21 shows the map of the 
bay area. The suitable tracing parameters were established to be the dissolved 
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oxygen, the electric conductivity , the turbidity and water color or 
transparency* A meter is being devised for quick quantitative measurement 
of the last parameter in situ. Figure 22 and 23 show the results for those 

parameters at the Cloquet intake. The times of poor water quality are 
evident from the graphs. To interpret these events, ve consider how the 
water quality at the intake could vary, by examining the Table in Figure 2h, 
The natural storm activity should increase turbidity, conductivity and 
color, but not the dissolved oxygen since lake water is saturated with dis- 
solved oxygen. Presence of the effluents in the lake water should on the 
other hand be characterized by an increase in conductivity and color and a 
decrease in dissolved oxygen leaving the turbidity relatively unchanged since 
effluents are relatively low in turbidity. Figure 22 and 23 show both types 
of events. It is seen that the latter type of event is displayed for the 
period of durable ice cover on the lake in late February. Odor in the lake 
water was detected at the same time. Pronounced odor characterizes the bay 
water. Lake water at other intakes was odor free* Temperature profile of 
the lake ruled out turnover. It is thus tentatively proposed, based on data 
shown in Figures 22 and 23, that cause of odor in late February 1973 at the 
Cloquet intake was due to St. Louis River effluents. Chemical tracing for 
next year should establish this better. • 

So far ERTS images for spring break-up are not available. It is hoped 
that subtle changes in images due to water color caused by the accumulated 
effluents can be discerned from the remote sensing data taken on days when 
the ice cover just begins to move out of the lake tip area. Should this be 
possible, an estimate of the extent of effluent accumulation could be made. 


i 
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Such information would establish the location of crucial sampling stations. 

I 

The remote sensing data supplemented by sampling measurements these points 
vill allow for calculation of the transport coefficient for effluents in 
the lake. Extensive sampling could also do the Job, however it would be 
prohibitive because of the large area that would need to be covered if the key 
sampling stations could not be predetermined. This year the safe ice cover 
lasted too short a time to collect sufficient data on sampling basis to 
analyze the wintertime concentration of effluents , at points other than 


the Cloquet intake. 
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92 v mhos/cm 

125 P mhos/cm 

D.O. 


2.5 p.p.m. 

' 0 p.p.m. 


Table Fig. 20. 
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Table II 


Normal Values for Water Quality Parameters 



Clean Water Values 
at Cloouet Intake 

Winter Values at 
Superior Bay 

Turbidity 

2 

6 JTU 

Specific Conductance (at 1.5°C) 

55 

110 p mhos/cm 

Color 

3 

300 Color units 

Dissolved Oxygen 

11 ppm 

2 ppm 


Calibration of Conductivity Probe VSI 33 S.C.T. 


KC1 Solution 

Temp. 

Probe 

Standard 

.0005 m 

23°C 

75 P mhos /cm 

73.9 P mhos/cm 

• 001 m 

2^°C 

1^0 p mhos/cm 

, 1^5 P mhos/cm 

• 01 m 

24°C 

1370 p mhos/cm 

ll 08 p mhos/cm 

\. 

Lake Water Station #1 

(150 ft. depth) 

90 p mhos/cm 

at 2^°C (old sample) 

Recent Cloquet Intake 

sample 

95 P mhos/cm 

at 24°C 


5U p mhos/cm at 1,8°C 


F* W. Q. Lab Intake 


3 ^ 


Ice Growth for Superior Bay • • 

1) Formation: The problem of formation of first ice may be separated 

into two periods. The first is the cooling of the body of vater 

to an isothermal state. The second is the cooling of the surface 
vater from this isothermal condition h°C ) to 0°C, and formation 
of a solid ice cover of ^ H in. The figure h in. was chosen be- 
cause once that thickness is reached, one is fairly assured that 
the cover is permanent. Experimental evidence shows that once the 
surface reaches 0°, the first h in. of ice form very rapidly. For 
the first stage of cooling, Newton* s Law of Cooling may be applied.* 

where T v « H^O surface temp. 

T a = air temp. 

k = Constant of proportionality 
t - time 

^ j 

let T w (t n ) = T n 
T w (t 0 ) = T 0 

* M. A. Biiello, Ice Prediction Curves for Lake and River Locations in Canada. 

U, S. Army Command, Hanover, New Hampshire. July 196^ 

B. Rodhe, On the Relation Between Air Temperature and Ice Formation in the 
Baltic, Georgrafiska annaler, vol. 3^, p. 175-202. 


dT 

dt ” 


dT v + kTydt = kT a dt 


e^aTy + ke kt T w dt * ke kt T a dt 


.ktr 

L w au = 

kt \ m at ^kt ' 


or d(T v e kt ) = T a d(e Kti ) 
Integrating from t = t Q to t 


T v e 


kt 


t~t T-l - 

t = t 


:n 


n = /“ T a d(e kt ) 
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: n e kt n - T 0 e kt ° = j T d(e kt ) 


T n ” + e“ ktn /* T a a(e kt ) 

t o 

Write integral as a sum by breaking (t 0 ,t n ) into n* equal short 
intervals , At ~ 1 day 


n 


T n = + e ~ kt “ r . v 

V=1 


T v (e ktlr -e kt v-l) 


vhere T y = Ave. value of T a on t v _-^ to t y , 
T "^V 

i.e. , T v = j- / T a (t)dt (At = Vtv-l) 


tv_l 


T n = T 0 e -k(t n- t o) + ” „ [ er k(t n -tv)_ e -k(tn- t v+At) j 

v=l ' 


t = T 0 e _k(tn-to) + (l_e~ kAt ) E Te' i(tn4v) 

V=1 

Selecting the period large, t n -t Q »1, 

r 0 e~ k ^ tn_t °^ = 0 . “ — 

To obtain a f omul a of iteration, replace t 0 by t n _ 2 * 

T n = + (l-e" kAt )T n 

or t£ = T'_ x + (l-^KV^) 
is an "adjusted air temp." 

Nov the value of (l_e~ kAt ) vas found by trial and error such that 
for n = no. days from starting point for to reach H°. 

The starting pt. vas July 1 (any day several months before freeze- 
over vould do) and vas taken to be the average temperature for the 
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month of June. This vas done on computer and found that 
(l-e - ^^) « .99 for the Superior bay. (This figure has 
not yet been calculated for the lake. ) It was assumed here 
that this atmospheric heat exchange was not affected by 
winds and currents. When the bay has become isothermal, 
the water surface temperature will be much more sensitive 
to air temperature. It has been found experimentally that 
the formation of solid ice cover occurs after a set number 
of accumulated degree freezing days, (32°F base temp.) be- 
ginning from when the body is isothermal. + This experimental . 
fact may be rationalized by assuming that during this second 
cooling stage the lake-water surface temperature is close to 
32°F, and the rate of heat loss is proportional to (T a -32)°F, 
where T & = air temperature (°F). This figure was 130 accumulated 
degree freezing days for the bay (ihxluth-Superior). Using the 
figures (l-e“^^) and 130 ADFD calculated from 1972 data, the 
formation of a solid ice cover was predicted within three days 
for the years 1969.5 1970, and 1971. The temperatures used in 
these predictions were the actual temperatures recorded, but in 
practice one would use actual temperatures up to the day the forcast 
is made, and then use forcasted and average temperatures after that. 
Thus one would be able to increase the accuracy of the forcast as 
each daily average temperature is recorded. This should allow for- 
casting freeze-over (to within five days) about two months (or more 
in advance. 

» 

t G. P. Williams, "Correlating Freeze-up and Break-up With Weather 

Conditions" Canadian Geotechnical Journal. Toronto V. 2, #4 1965 
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2) Growth: There are two major types of ice growth, black ice, which 

is formed from the freezing of lake water, and white ice, which is 
formed from the freezing of slush on top of the black ice* The a- 
mount of black ice (as a function of time) may be shown to be pro- 
portional to the square root of ADFD (accumulated degree freezing 
days)* Consider a unit area of ice, thickness x. 


T„ 


j£L 


T 


* ^ dx 

Tf / 




T: 


— « (t-T ) 

dt “ v 

Let ~ * (T.-T) = k 2 (T-T ) 

dt x 1 ^ a 

where k^, kg are constants. 


X dx 1 dx 


at + dt = (T i" T) + (T " T a) = V T a 


j- dx + ih a* = (T.-Tjat 


...kg 


x=s 

L + fe) = / . < T i- T a )dt 


1 

t«t" 

/ 

t=t y 


‘ \ 


where x(t > ) = S (the ice thickness is S at t^t^) and t is taken as 


the first day where ADFD(t) stays positive. 

3^ + % “ j ( V T a > d * = MFD 


For forced convection, « k g (due to motion of air over ice) 

t' • - 

S 2 - 2k ■, / (T.-T )dt or S « [ADFDj 1 ^ 2 

t o 

1 /2 

For Superior Bay, x(t) = (5.6")[ADFD3 . (5.6" determined 

experimentally) 
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Since wintertime acquisition of data required ice penetration, the 
collection and examination of ice data itself could hardly be avoided, both 
for curiosity and safety reasons. Excellent data for one point on the bay 

ll tc 

vas available at the U. S. Army Corps of Engineers. Since our early 
measurements or ice cover indicated that 1972-73 ice growth vas unusual, 
it vas interesting to attempt to predict the freeze-up times and the break- 
up time, for this year using the available data. This year for moderate ice 
thickness (0-20" ) the grovth of black ice vas veil correlated with the ac- 
cumulated degree freezing days. The ice thickness vas proportional to the 
square root of degree accumulated degree freezing days as can be seen from 

Figure 25. The reason the formula holds so veil came about from the fact that 

6 - 

this year negligible white ice had formed* The departure from this relation 
past 20" in thickness could be attributed to river currents which are not 
steady on account of the lake, and the presence of the two entries.. Also 
the effects of the lake ice cover on the air temperature on the bay and the 
unusually mild weather itself could account for^this departure* It is seen 
from Figure 2 6 that the relation is not valid for years when the accumulation 
of white ice vas not negligible. The prediction of ice breakup would depend 
on the ice type since thawing curves like the ice growth depend on the nature 
of the ice. Black ice shows regular grovth and decay. Analytical description 
of formation of white ice and its effect on breakup time would be difficult at 
best. To predict the breakup time however, some effective description, 

U. Private communication, Clarence Wong, U. S. Army Corps of Engineers at Duluth. 

5. Mr. T. D. Brennan, U. S. Coast Guard. 

6. J. B. Shaw, McGill Sub Arctic Research Laboratory. 
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accounting for white ice, was necessary. Records of plot of total ice 
thickness as a function of time for the past four years indicated that 

the turnaround tine, for ice growth on the Superior bay occurred on or be- 
fore March 6 . Thus a prediction of breakup based on ice thickness for that 
date could be made using average ice decay rates. The question arose as 
to what was the effective thickness of ice cover on that day, realizing 
that the properties of black and white ice so far as decay rates and 
mechanical strength are concerned, are quite different. As an approximation 
of an effective thickness, we took the thickness of black ice plus one half 
of the thickness of white ice. Often this value corresponds to ice thickness 
below water level, but the water level fluctuates on account of ice buckling. 
Upon assigning to white ice one half the effective thickness in terms of black 
ice, and averaging over successive readings for each year ve obtained, Figure 
27 which shows the total effective ice thickness for .the past five years. 

Once the turnaround in ice thickness occurs, it is seen that the. ice 
decays at an average rate of .227 inches per day until it reaches a critical 
value. At such time parts of the bay where currents are large, open up, while 
the ice cover on Superior bay channel honeycombs and cracks, presenting little 
obstacle to ship traffic. This critical ice thickness for the mid area of the 
Superior bay is 15”. Thus if one projects the decay of ice from the March 6 
thickness, the latest maximum thickness date, to the 15” thickness , one ar- 
rives at predicted breakup as shown in Fig. 28. Since for 1973 the maximum 
in the ice thickness occurred more than a month early, it was interesting to 
use this method to predict the harbor opening of this year. Fig. 28 shows 
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that the predicted date was March 22, which was correct, though shipping 
did not start until March 29. The shipping however depends on other factors, 
such as opening of the locks and the ice packing at the tip of the lake. 

The measurements on the ice cover at the tip of the lake were made 
during the times of collection of water samples on the lake at points inter- 
mediate to the Cloquet intake. The ice cover on the lake was, however, 
minimal this year. It was hazardous for all except two weeks time. The ice 

T 

.cover on the lake was treated according to procedure of T. L. Richards. 

Since the tip of the lake is open to the rest of Lake Superior, it could not 
he treated altogether as a separate "body of water. If no major storms were 
present, the tip area could act like a shallow lake* Using characteristic 
parameters of antecedent heat, and accumulated degree freezing days, and 
the assumption that the tip of the lake behaved like a shallow lake, it 
is estimated that 100$' cover would develop after about 800 accumulated 
degree, freezing days, which would yield (averaged over records for 1967- 
73) the' 100$ freeze up date as January 15 and agrees well with projected 
growth of initial cover. However, due to wind storms, partial ice cover 
at the tip of the lake is often removed and new warmer water is brought 
into the area, initiating in effect a new ice growth period, Rig. 29# 

After an accumulation of 1000 degree freezing days , the entire lake 
cools sufficiently so that this factor is less important. Thus a durable 
cover over the tip of the lake forms after about 1,100 accumulated degree 
freezing days which falls at the end of January. Fig. 29 shows the ice 
cover for 1972-73,' using ERT3 data for the winter. The 100$ area considered. 


7. T. L. Richards, Monthly Weather Rev.', Vol. 92, 6, 297- 



is the tip of the lake from Knife River to Duluth* The^ ice thickness after 
freeze up seems to display growth proportional to the square root of the 
degree accumulated freezing days. Ice thickness data for this part of the 
lake is unavailable for other years. The main problem to navigation in so 
far as. the ice at the tip of the lake is concerned, comes from ice packing 
due to ice driven into the tip area in the early spring by N.E. winds. This 
phenomenon for instance caused considerable difficulty to shipping in 1972, 
when the ice was piled up to 10 ft. depths, and remained in the area until 
June 9» 1972. Similar event caused l8 M pack to develop late in March of 
this year. The extent of packing depends on the ice cover on the entire 
lake. Further consideration of this problem is planned with aid of the 


ERTS satellite data. 
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Our measurements on the effects of the experimental bubbler vere per- 
formed as supplementary information to the data collected by the National 
Biocentric Corporation » which vas charged with assessment of its operation 

for the U. S. Army Corps of Engineers. Five sampling points vere established 
in the area of bubbler as shovn in Fig. 30. The measurements of vater quality 
parameters at these points vas used^for the evaluation purposes. Intensive 
measurements vere made for three sampling days, one each before, during and 
after the operation of the system. Ve also took data on several other dates 
as part of our investigation of the bay* The auxiliary stations used for 
our ovn purposes are shovn in Fig. 30. 

The assessment of the operation of the bubbler so far as ice cover was 
concerned, vas based on the ice growth* data for points adjacent to the bub- 
bler. 

Fig. 31, curve 1, shows the ice growth for area upstream from the bubbler. 
Curve 2 on the figure shows the estimated growth of ice over the bubbler area 
if the area was not disturbed. This vas done by using the same functional de- 
pendence for ice growth over the bubbler area as for station #5* Initial values 
for ice thickness were based on measurements made before the bubbler was in- 
stalled. Curve 3 shows the actual ice values over the bubbler affected area. 

Two discontinuities occurred; one on January k because of a wind storm vhich 
caused considerable surge of lake water into the bay, breaking up the ice 
cover; and one on January 15 because of unusually mild weather, at that time 
the bubbler perforated the ice, creating open water right over the bubbler 
outline. For more* normal conditions the bubbler area, including the section 
right above the bubbler outline, vere frozen over. The ice right above the 
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bubbler outline vas, however, substantially weaker and remained in that 
state even after the bubbler, operation ceases. 

When the bubbler vas shut off, the ice cover quickly returned to near 
normal projected thickness. The veaker area directly over, the outline of 
the bubbler formed a pressure ridge about 2 ft. high. When the thawing 
period began this ridge area deteriorated first, shoving open water a veek 
earlier than adjacent areas. 

Based on this years performance then, (year unusual in character) it 
is estimated that the breakup of ice in the bubbler area vas speeded up by 
one to two veeks on account of the weakened section formed directly above 
the bubbler. Even in the absence of the weakened ice and the subsequent 
pressure ridge, one could estimate that the ice cover over the bubbler should 
deteriorate sooner. It displayed many imperfections in growth, (characterized 
by trapped air and discolorations giving appearance of entrapped algae growth) 
and its final thickness was about 2 inches lower than projected normal thick- 
ness. This according to the thawing curves. Fig. 28, would make the bubbler 
area navigable 8 days sooner than the adjacent area on. the Superior bay. 

The assessment of the effects of the bubbler on water quality was dif- 
ficult. The location of the bubbler vas chosen for its anticipated aid to 
shipping as such the choice was excellent. The area does, however, have the 
problem in that it is not a homogeneous body of water. Fig. 32 shows that 
the bubbler is located at the junction of the muddy Namadji River and the - 
much less turbid St. Louis River. The proximity to the lake entry also com- 
plicates the problem as the lake water periodically flows into the bay. Some 
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assessment of the effects of the bubbler on water quality can be made, 
however, particularly if comparisons are made for selected sampling stations. 

We limited most of our measurements only to those physical parameters 
which could be easily taken in situ, and for which ve felt pertinent relations 
to water mixing could be established and discussed* (Other measurements were 
performed by Biocentrics, Inc.) 

Basically, one would like to assess the effect of the bubbler in prevent- 
ing the natural settling of pollutants, and on its possible reintroductian 
of settled out material into the water. The question in this case was 
whether the mixing was largely due to the bubbler or due to the natural 
mixing processes at the junctions of the rivers and the lake entry. 

Aerial data was most useful in this "case. Fig. 32 gives good indication 
that the site #1 on the bubbler should reflect the properties of the St. Louis 
River. One would expect the same from site #5* These two sites could then be 

compared. The station at the tip of the bubbler, off the old ore dock, could 

v % 

serve as an indicator of lake water surging into-the bay, the surges do indeed 
occur, as evidenced by the water current observations at the Superior entry, 
and water percolating through the ice cracks in the slip of the old ore dock. 
It is anticipated that the prevailing water flow in the area; except for areas 
right next to physical protrusion, such as the U. S. Army Corps pier; is 
normally laminar. This is somewhat borne out by the results obtained for 
the times of sampling after the cessation of bubbler operation. 

The data in Figs. 33-39 show that a normal profile with depth for the 
conductivity and the dissolved oxygen had developed by December 13 for the 
Superior bays. Considerable settling is shown by the depth profiles for 
the site off Barkers Island, which is downstream from the Superior sewage 
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plant. The first intensive sampling day was December 21, Fig. 3**, 35 • 

The bubbler area vas open. The tugboat, Illinois, courtesy of Great Ijakes 
Shipping Company, vas used in collecting the data. The data, for December 

21 shows that the water in the bubbler area is veil mixed, particularly if 
comparison of the shape of conductivity profiles is made for sites #1 and 
#5* Furthermore, the average conductivity of water was substantially higher 

at site '#1 indicating extra conductive suspended material either due to 
bubbler action or the vessel traffic. Site #5 was not appreciably dis- 
turbed by the tugboat on the sampling day* The results for subsequent days . 
show water quality parameters for Barkers Island site and the upper channel. 
The upper channel displays little variation of conductivity with depth. Gen- 
erally pronounced conductivity profiles were present primarily in the Superior 
bay and the St. Louis Bay where major sewage and industrial waste sites are 
found. The next day of extensive measurements vas January 15 » Fig. 36, 37- 
The data for the day shows again considerable mixing at site #1 as compared 
to site # 5- The ship traffic had stopped by January 1973, and water 
current was downstream on that day. The day vas unusually warm and the ice 
directly above the bubbler vas either well perforated or vas shoving open 
spots. 

The oxygen over the bubbler area had increased ^ 2 ppm on the average. 
This would depend mostly on the rate of flow of water over the bubbler. The 
oxygen profile vas determined using proble measurements and Winkler (Alkali- 
Iodide Ozide) analysis on samples every 5 ft. depth for the extensive measure- 
ment days. The comparison of the conductivity curve for site #1 and 8 5 shows 
that the average conductivities for the sites are nearly the same, 
unlike the case for December 21. 

ft 
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One could then deduce that the water at site 85 "become well mixed by the 
time it reached site #1, but if the increase in average conductivity on 
December 21 was due to bottom sediments being resuspended , then by January 
15 the process was no longer tailing place and though the bubbler was pre- 
venting settling, it was not adding new conductive bottom material to the 
water on that date. The station off the tip of the bubbler reflects the 
parameters similar to those at station 85 while station 83 seems to reflect 
the regular winter conditions at that point which changed little later on. 

In the areas away from the bubbler, a slight temperature gradient had 
developed by that date. This was typical of the temperature profile later 
on, and may veil be of considerable importance in improving the effectiveness 
of a bubbler operating in mid-winter. The profiles for January 30 and 
February 6 displays depth variation in the parameters after cessation of 
bubbler operation. Some natural mixing is indicated on February 6. The 
data for February 13 shows the results for the times well after the '..operation 
of the bubbler has stopped, so does the data’ for later, days. -It is noticed 
that for times after cessation of operation the bubbler area generally dis- 
played stratification indicative of the upstream conditions on the Superior 
bay. Station 8 1 and #5 show small differences , particularly in the tempera- 
ture profile. Natural mixing may have been the cause of this* The effect 
is relatively small however compared to mixing indicated during the bubbler 
operation. The available data thus suggests that the bubbler introduced 
1-2 ppm dissolved oxygen into the water at the site ft 1 and sufficiently mixed 
the water to distribute . the conductive ions almost uniformly with depth. The 
initial operation monitored on December 21 seems to have also caused 
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substantial increase in the average conductivity, of the vater,. perhaps by 
resuspending the settled out particles. However, the cause of this cannot 
be established since considerable turbulence was caused by vessels operating 
in the area. A more intensive study for a bubbler, located further away from 
the entry, should be made particularly with regard to the effects of bubbler 
on the bottom material. “The operation should also be extended to the entire 
vinter to sea how the bubbler operation improves with the natural development 
of temperature gradient in the channel. This could shift the ice breakup to 


a much earlier time. 
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Summary of progress : 


llio studies for the past nine months have established the type of 
measurements vhich should be pursued in subsequent investigation of the 
year-round vater transport in the western tip of Lake Superior and the 
St, Louis River bays. 

Following is' our plan for future studies. 

"V 

Objectives : 

Summer-Fall 

1. Continuation of studies of the upvelling along the N.E. shore 
and the eddy currents at the tip of the lake with aid of ERTS 
images shoving the turbidity plumes and patterns at the tip of 
the lake. 

2. Continued investigation of vater quality in the bays and in- 
vestigation of the transport of effluents out of the harbor • 
area into the lake and hov these pertain to the vater quality 
at the tip of the lake* 

3. Investigation, of the mixing processes in the Superior bay - 
using the established parameters suitable in detecting 
mixing and the vater current measurements at the entries. 

A new instrument for current measurements vill be set up 
by the U. S. Army Corps of Engineers at both entries. 

Winter 

1. Study of vinter currents in the lake through physical and 
chemical tracing of the effluents. 
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2.. Study of ice growth, cover and ice packing in the lake tip 

* area using ERTS images; end continuation of harbor ice 

studies together vith the environmental effect of bubbler, 
should these be in operation again. 

Plan of operation : 

Summer-Fall 

To establish further the nature of the upvelling along the North Shore 
for westerly winds, horizontal temperature gradient, and accompanying water 
velocity profile with depth will be monitored at Knife River site and a site 
in the vicinity of the Cloquet water intake. 

To obtain quantitative data on the eddy current depth profiles for 
temperature, turbidity, water velocity, and the thermocline activities will 
be made at a site near the Cloquet intake and a point 6 miles directly N.E. 
from the station where the ERTS images show a northerly turnaround in the 
currents. The data will be then correlated vith the ERTS data on turbidity 
plumes . 

To produce quantitative data sufficient for analytical study of the 
eddy current, at least five continuously monitoring stations should be 
established at the tip of the lake. If our application to other agencies 
for the support of more extensive ground truth data collection is accepted, 
such station will be established. Under present conditions of manpower and 
instrumentation, however, only two regularly visited sites can be maintained. 

During the times of rough sea conditions on the lake, the work on 
monitoring of water quality in the Superior bay and the St. Louis bay will 
be continued on the same basis as in the past. The data will be examined 
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for influence of effluents on the vat er quality at the Cloquet intake, 
which vill he continuously monitored by the city of Cloquet. 

Two intensive 3 day studies consisting of continuous monitoring of 
vater quality in the Superior bay and the vater currents at the entries 
vill be made to evaluate the mixing process in the bay. This material 
together vith other measurements in the bays vill serve the purposes of 
masters thesis for Mr. Kirby Stortz. A detailed operation plan is a- 
vailable should NASA require it. 

Winter 

Winter studies vill be maintained at the 1972-73 level. The chemist ry 
department at UMD vill have complete their operational setup of the experi- 
ment this year so that the tracing of effluents by chemical means vill add 
to our information of wintertime vater transport in the lake and vill identify 
further the cause of odor at the Cloquet intake. The local U. S. Army Corps 
of Engineers in going to install a year-round system for vater current 
measurements at the entries. This should substantially aid us in studies 
of mixing processes in Superior bay. We had requested support for our 
wintertime studies from the EPA as we felt that the studies were germane 
to their interests , since three municipal vater intakes are located at the 
tip of the lake* The agency has rejected our request. 

Studies of the onset of ice cover and its growth in the western arm of 
Lake Superior will be pursued further. ERTS data for the ice cover vill be 
used in relating the ice cover to the spring packing of ice at the tip of 


the lake. 
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As a part of our understanding of KASA’s practice to serve the 
community , we enumerate the possible application' of our work for the 
benefit of the community: 

1. Work on water transport in the winter and summer is necessary 
to water users of the area. Currently all intakes experience 
some difficulty with water quality, though it is intolerable 
only at the Cloquet intake. That intake must either be re- 
located or a filtration must be provided designed with sufficient 
knowledge of water quality parameters for any further degradation 
in water quality at their intake. Such could for instance arise, 
should extensive bubbler systems be installed in the Superior 
bay for purposes of extension of shipping season. 

2. The knowledge of water transport is also essential for purposes 
of determination of the effects of dumping of harbor dredging 
in the lake. In the past dredgings were disposed of at the tip 

' of the lake. Such -practice will resume in the future when the on 

" ^ 

land disposal sites are unavailable or too costly. It is antici- 
pated that this will happen within the next 10 years because of 
limitations of the on land sites , and likely much sooner for 
economic reasons. By that time the ongoing program in pollution 
control should remove the dangerous pollutants from the harbor. 
However, there will always be considerable amounts of organic, 
fibrous material which could degrade the water at the intake if 
proper dumping sites are not selected. Dr. Keith Larson of the 
U. S. Army Corps of Engineers at Minneapolis is well aware of 
this problem and is considering our work as a part of the Corps 
long range program. 
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The winter studies on the bubbler operation and the work on 
ice conditions in the 'bay and the lake are useful, to the 
shippers. Excellent program, for monitoring overall ice cover 
on the lake exists and will improve with ERTS program. How- 
ever, quantitative data on ice growth at the tip of the lake 
is sparse. The study of ice conditions in the harbor and, ice 
growth and packing at the tip of the lake is important to the 
proposed extension of the shipping season. 

Involvement in community : 

1. We have had close cooperation from the Cloquet City Engineer, 

Mr. Bruce Boyer, in use of their facility and exchange of in- 
formation on their problem with water quality at their intake. 

2. This year we will be involved to a limited extent with the 

harbor dredging program through the National Biocentrics 
Corporation and the U. S. Army Corps of Engineers. . ■ 

3. Our winter studies on bubbler have been’ inade in cooperation 
with the Biocentric Corporation and the local U. S. Army Corps 
of Engineers, the local office under Mr. Clarence Wang. 

The monitoring of water quality on St. Louis River and its bays 
is a cooperative program with Mr. A. Biele, the chemist at the 
Lakewood Pumping Station for Duluth. He is responsible for 
most of the chemical analysis of the sample's. 

I wish to thank Mr. Clarence Vang; Mr. Bruce Boyer, and Mr. Art Biele 
and their personnel, for their advice and cooperation on this project. I 



vould also like to thank Dr. Ted Odlaug and Dr. Ted Olson Tor their 
discussions, guidance, encouragement, and the use of their research 
facilities. Finally, I vouid also like to thank my students, Kirby 
Stortz and James Rohlf for helping me out in the task of data collection. 
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Artificial defoliation in corn was used to explore the 
useful ucs* of ;u*ri:d photography in deter ting crop-insert 
infestations. Defoliation on top of plant was easily detected, 
while that on the base was loss so. Aero infrared film with 


W ratten SOB filter gave the best- results, and morning 
flights at the. scale of 1 : 1584(1 arc recommended. How 
direction, plant growth stage, and time elapse sinre 
defoliation were not important factors. 


Karvae of both army worms, Pseudalctia nnipuncla 
(Haworth), ami grasshoppers defoliate rorii plants. 
Information on the amount of leaves destroyed is useful 
in assessing potential crop losses, and assists in es- 
timating the levels of pest populations involved. To 
obtain such information by ground checking requires 
considerable labor and time. Time is particularly crucial 
because early availability of this information is essential 
for preparing plant-protection measures. 

Aerial photography has been used in various aspects 
of agriculture (Brcnchley 19G8). Work has been done on 
forest insects (Heller et al. 1959, Baj'zak 1966, Wert and 
Roettgcring 1968, Harris 1971); range insects (Meyer 
and Woolfolk 1967); ami fruit insects (Hart and Myers 
1908). No information is available on detecting ccreal- 
erop-insect infestations. 

We explored the usefulness of aerial photography 
in this latter respect. The work was done on the Rose- 
mount. Experiment Station of the University of Minnc- 
sntain Rosomount in 1970 and 1971. 

General Methods.— A erial Photography and Image 
Analysis . — Multibnnd reconnaissance equipment of the 
University of 'Minnesota, College of Forestry, in the 
form of a 70-mm quadrioamcra unit and a photo aircraft, 
were used. The film-filter combinations, photo scales, 
overflight times, and film enhancement procedures used 
in 1970 ami 1971 follow. * 

1970 1971 


Film/fiUer combination 


Ektachrome MS/Wrattcn 2A 


X 

Ektachrome Infrared/Wratten 12 

X 

X 

Ekt-acolor/Wratten 2A 


X 

Panchromatic Plus-X/Wrattcn 8 


X 

Panchromatic Plus-X/Wratten 58 

X 


Panchromatic Plus-X/Wratten 25A 

X 


Aero Infrared/Wratten 89B 

X 

X 

Photo scale 

1:4000 

X 


1 :6336 


X 

1:12000 

X 


1:1.5840 


X 

1:31680 


X 

Time of flight (OST) 

0900-1000 

X 

X 

n oo-i 200 

X 


1300-1400 

X 



1 Pi»p*rr no. 8070 flrfentific Journal Series, Minnesota Agricultural 
Eiprriment Station, Bt. Paul, Minn. Received for publication Aux. 18, 
1072. 

» Profcwor, Department of Entomology, Fisherle*. and Wildlife; 
formerly Research Fellow; and Profeeaor, College of Forestry, Ifni* 
vemity of Minnesota, respectively. 


Film enhancement procedures 

PS Mini-Addcol X X 

PS Digicol Viewer X 

Photo-optical Density Slicing X 

ISl VP-8 Image Analyzer X 


To insure that all possible photographic characteris- 
tics were considered, both mnltiband ami Ektachrome 
infrared transparencies were subjected to various types 
of enhancement (Meyer and Chiang 1971): (1) the 
1970 70-mm multispcctral (three) positive transparen- 
cies were optically combined through component filters 
by means of the engineering prototype of the PS 
(International Imaging System) Additive Color 
(Addcol) viewer; the current model 6030 Mini-Addcol 
combiner was used for analysis of the 197 1 multispcctral 
imagery; (2) the 1970 Ektachrome infrared trans- 
parencies were enhanced through a process of color 
separation using a photo-optical density-slicing tech- 
nique; the resulting isodensity image (10 levels) was 
then reproduced in 3M Color Key, registered, and 
mounted; (3) one of the better Plus-X/Wratten 25A 
1970 season transparencies was evaluated with the I 1 2 S 
4000 Digicol viewer; (4) the Jul. 31, 1971, Aero In- 
frared/Wratten 8911 transparency was enhanced on the 
ISI VP-8 Image Analyzer. 

Field Plots . — Two fields, P-10 and A-15, were used 
in 1970, and one, P-10, was used in 1971 (Fig. 1). Tho 
general area has Waukegon Bill loam type of soil. Some 
characteristics of the crop and planting, and information 
on defoliation follow. 


Row direction in the field 

1970 

1971 

E-W (field P-10) 

X 

X 

N t S (field A-15) 

X 


Row ividth 

Plant population (no. plants 
per acre) 

30 in. 

30 in. 

Field P-10 

19,250 ’ 

20,666 

Field A-l 5 

21,000 

Com variety 


Field P-10 

Pfistcr 

Pioneer 


SX-3935 

3784 

Field A-15 

Funks 


Defoliation levels 

(a) 0 leaves removed 

4222 


(control) 

X 

X 

(b) 2 basal leaves removed 


X 

(c) 4 basal leaves removed 

X 

X 

(cl) 6 basal leaves removed 


X 


Je) 8 basal leaves removed X 

.0 2 top leaves removed X 
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(g) 4 top leaves removed 


X 

(h) <» top leaves removed 

(i) 4 basal leaves removed 


X 

on alternate plants 


X 

()) 4 alternate leaves 



removed 

Numl>er of replicates in each 


X 

field 

6 

2 

Plot size 

24 rows X 

24 rows X 


75 ft. 

75 ft. 

Subplot, size 

8 rows X 

8 rows X 

Dates of defoliation/no. ex- 

25 ft. 

25 ft. 


tended leaves at that 

time Aug. 14/16 July 10/11 

Aug. 21/18 Julv 20/13 
Aug. 3/16 

Dates of flights/days after 

defoliation Aug, 24/13 July 20/1 

Aug, 24/10 July 27/1 
and 8 
July 31/5 
and 12 
Aug. 4/1, 

9 and 16 

PitncKPUHKK and Results. — 1970 Study . — Artificial 
defoliation was to be done late in July, the usual time 
when army worm larvae invade corn fields under 
Minnesota conditions, but it was delayed until Aug. 14 
and 21 because of weather problems and aircraft avail- 
ability. Flights were made Aug. 24. The plants had 
reached their maximum height with approximately 16 
leaves, ami the silk was in the drying stage. Two levels 
of defoliation were made* 4 and 8 basal leaves were 
removed, both simulating armyworm larval feeding. 
The 2 levels and a control were each replicated 3 times 
within a plot. All leaves removed from the plants were 
hauled off the plots. Flights were made in the forenoon, 
at noon, and in the early afternoon. 
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F 10 . 1. — Fields in 1970 and 1971, showing types of 
defoliation (a-j, see text). 


The transparencies were examined both with unaided 
vision and under magnificat ion. Du all film-filter t om- 
binations, and at both scales, dejoliation level e was 
dearly visible with or without magnification fi'ig. 1). 
This iiigh degree of visibility was duo to 2 hietors: (a) 
stress caused by removal of approximately \-± of the 
leaves on each plant, ami (b) somewhat increased 
vertical exposure of the soil. In no case was defoliation 
level c visible on the original photographs. This was 
not particularly surprising in view of the probable lack 
of effect of this relatively light defoliation upon plants 
so near maturity. 

The 3 enhancement procedures described earlier were 
applied to the transparencies. None made defoliation 
level c visible. However, with the 1*S Add col viewer, 
the visibility of defoliation level c was improved, and 
other features such as crop maturity differentials, 
drainage pallerns, compaction due to past animal 
trailing, were realized. With the H.S 4000 Digirol 
viewer, details discernible on the. original imagery 
became much more visible to the unaided eye. 

The results showed several points of methodological 
importance: (1) row direction, E-W or N-ri, did not 
influence over-all pictoral quality and rendition; (2) 
the quality of forenoon photography was the best, 
afternoon the next, and noon the poorest; (3) the 8-row 
plots were large enough for experimentation; (4) 
replications within a block were not necessary. 

1971 Study . — Since 1970 results showed that E-W 
and N-S row directions did not make any difference, 
only 1 field (P-10) was used in 1971. We photographed 
only in the forenoon and replications were reduced to 2. 
Hut more types of defoliation below the level of 8 leaves 
removed were studied: 2, 4, and 6 basal leaves were 
removed, simulating armyworm larval feeding; 2, 4, 
and 6 top leaves were removed, simulating grasshopper 
feeding; and 4 basal leaves were removed on alternate 
plants, and every other leaf on the plant (4 leaves in 
total) was removed. The last 2 types are to simulate 
lighter levels of insect infestation. 

Defoliations were done July 19, July 20, and Aug. 3, 
when plants had 11, 13, and 16 leaves, respectively. 
Tassels were out Aug. 3. Flights were made July 20, 
27, 30, and Aug. 4. Thus photographs were taken at 3 
stages of plant maturity, and from 1 to 16 days after 
defoliation. 

Transparencies were viewed under stereomagnifica- 
tion. As shown in Fig. 1, there were 9 subplots with S 
different types of defoliation and a control. The subplots 
showing colorations lighter than the field in general 
were noted. The lightest was assigned a rating of 1, the 
next lightest 2, and so on. A subplot did not have a 
rating if it was not distinguishable from the field in 
general. In a given plot, there might be 2 subplots 
showing the same degree of lightness, and they were 
assigned the same rating. Thus in a plot there could be 
more than one subplot rated 1, or 2, or 3. It should be 
noted that the subplot rated 1 in a given plot might not 
be of the same lightness as that in another plot. 

During the process of evaluation, the levels of de- 
foliation were not made known. Two persons viewed 
the transparencies and made evaluations independently, 
and produced similar results. 

Tabic 1 includes data to compare the 4 film-filter 
combinations and three time lapses at the photo scale 
of 1:15840. The film-filter combination Aero Infrared/ 
Wratten 8911 gave the best quality. More subplots were 
distinguished (see figures in the right hand column) 
with that combination than with the others. There 
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T.Wc 1. — Ratings of different lypes of defoli.lion with 3 time lapses, and 4 liim/fllter combinations ( pboto scale 1 : 'SHIP) 


Dnlc of 
defoliation 


Data of No. days 
flight lapsed 


Aug. 3 Aug. 4 


July 26 Aug. 4 


July 19 Aug. 4 


Film/filtcr 
mm hi nation 


Arm Infrared/ 
Wrattan 8911 
Kktiichrume Infrared/ 
\Y rattan 12 
lOktarhronin MS/ 
Wrattan 2A 
Kktarnlnr/ 

Wrattan 2 A 
Aero Infrared/ 

IV ratten 8911 
Kktarlnome Infrared/ 
Wrattan 12 
Kktiichrnmc MS/ 

\V ratten 2A 
Kktacolnr/ 

Wrattan 2 A 
Aero Infrared/ 

\V ratten 8911 
Kktaehrome Infrared/ 
\V ratten 12 
Kktaehromc MS/ 

W ratten 2 A 
Rktaeolor/ 

Wrattan 2A 


were no consistent differences among the 3 time lapses 
of 1, 9, ami 16 days. 

Table 2 is organized to compare the effects of 3 
factors on the photographic resolution of the different 
levels of defoliation. The 3 factors arc: (I) degree of 
plant maturity as reflected by the date on which 
defoliation was made, (2) time lapse as reflected by the 
number of dnys between (iefoluitinn find flight, iUiu (^) 
photo scale, the data include only the film-filter com- 
bination of Aero Infrared/Wratten 8911, since this 
gives the best, photographic, quality. Unfortunately,, 
mechanical problems during the scale III Might on 
July 20 ruined the photography. 

Comparisons can most readily be made by examining 
the number of times the subplots received ratings. The 
figures in the 1st and 3rd columns on the left reflect the 
effect of plant growth stage when defoliation was made 
(1st column) and the time lapse between defoliation and 
flight (3rd column). Figures in the last column on the 
right show no consistent differential result from these 

2 factors. 

The figures in the next to the bottom row reflect the 
effect of photo scale and type of defoliation. Scales 
1 :6330 and 1:15840 were about equal in clarity, while 
1 :31G80 was less clear. The bottom row gives the overall 
dotoctahilitv of the types of defoliation based on . fre- 
quency of being rated and values of the ratings. The 3 
types most consistently detected, in order of coloration 
contrast are h, g, and j. The next types arc d and f. 
The least detected were b, c, i, and a. The differential 
among these was slight and accuracy of detection very 
low; as reflected bv thr fact that level a, which is the 
control with no leaves removed, was given a rnting on 

3 occasions, although all at low levels of lightness. It 
should 1m* mentioned, however, that there was no such 
confusion at the scale of 1 :6336. 

PnACTiCAL Implications. —The 2 years’ results have 
the following implications which are of practical sigmn- 


Suhplots 

(defoliation levels) 

a b e d c f ghi.j 


No. subplots 
rated lighter 
than general field 


3 4 4 8 


1 1 
3 2 I 
2 2 1 
3 3 2 1 


cancc in using aerial photography for detecting de- 
foliation: 

1. Defoliation in corn may be detected with 
reliability regardless of row -direction. 

2. Detection is feasible within one day after 
defoliation; it is not necessary to wait for intensive 

plant stress. \ 

3. Detection is feasible when plants are at the 

11-leaf stage througKtasscling. , 

' 4. The scale of 1 :6336 gave most reliable detection, 
but considering the economics, the scale of 1:1 5840 is 
to be recommended. ^ 

5. Aero Infrared film with Wrattan 8911 filter 
gave the best results among the 7 combinations tested. 

6. Morning flights arc recommended. 

7. Defoliation on top of plant is easier to detect 
than that on the basal part of the plant. 

The artificial removal of corn leaves in square plots 
probably is easier to detect amid normally growing corn 
fields than natural infestations. It is probable that corn 
carworm or armyworm infestations would be more 
uniform over corn fields, ami the contrast levels would 
be more difficult to define. One other factor which was 
not involved in this test, because leaves were artificially 
removed, is the likelihood that grasshoppers or army- 
worms would not consume entire corn leaves. Partially 
consumed leaves would probably change color and Ik: 
better detected with normal color or color inf rami film 
than bv infrared film with an 89 B filter. 

Our next step will be to test the potential of the 
method in fields with real armyworm ami /or grass- 
hopper infestations. The information is published now 
with the hope that the method can be tested by en- 
tomologists elsewhere when the opportunity arises. 

Acknowlkikimknt. — Funding for the final image 
enhancement and analysis by NASA’s Office of Uni- 
versity Affaire (Proj. NGL-24-C05-263) is gratefully 
acknowledged. 



Table 2.— Ratings of dilTercnt levels of defoliation with 3 scales, 3 plant 


Scale I (1:6336) 
sub-plots (defoliation levels) 

Date of Date of Xo. days 

def. flight lapsed Rep. Hot abed fgh i j 


7/19 

7/20 

1 

I 

F 


3 


2 

1 

■ 

3 




II 

D 


2 

3 

2 

1 


2 

7/26 

7/27 

1 

I 

E 


3 


2 

1 


3 




II 

C 


2 


3 

1 


3 

, 8/3 

8/4 

1 

I 

B 


2 


3 

1 


3 




n 

A 


3 


2 

1 


3 

7/26 

7/31 

5 

i 

E 

3 

3 

3 

2 

1 

4 

3 


• 


ii 

C 


2 


2 

1 


3 

7/19 

7/27 

8 

i 

F 




2 

1 


3 




n 

D 


3 


2 

1 


3 

7/26 

8/4 

9 

i 

E 


3 

3 

2 

1 


3 


. 


ii 

C 


3 


2 

1 


3 

7/19 , 

7/31 

12 

i 

F 


4 

4 

2 

1 


3 




ii 

D 


3 

4 

2 

-1 


3 

7/19 

8/4' 

16 

i 

F 




2 

1 


3 




ii 

D 




2 

1 


3 

No. sub-plots rated lighter 

than general ficld b 


0 0 1 

13 

5 16 

16 

1 

16 

Overall rating 





9 8 7 

4 

5 

2 

1 

6 

3 


* Only Scales I and It are included. Others include Scutes I, II, nnd III. 
*> Scale* I and It include S groups, Scale III includes 7 group*. 
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June 1973 Chiano et al,: Aerial Photography op Defoliation in Corn 
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Fio. 2 . — Plate A, Optical combiner view of the 1970 Plus-X/W ratten 58, Plus-X/Wratten 25A 
Wratten 89B transparencies of Plots I Oeft) and II. Plate B, Digicol Analyser view of the portion of the Plus-X/Wratton 
25A transparency showing Plot II. Original negative scale ™ 1 ! 4,000; north is at top of page. 


i 



Fio 3.— Plate A. Optical combiner view of the July 31, 1971, Plus-X /Wratten 8 and Aero Infrared/Wratten 89B 
transparencies of Plots C. D, E and F (from left to right). Plot F is only partially shown. Plate B, VP-8 Image Ana- 
lyzer view of the Aero Infrared/Wratten 89B transparency of Plots C, D, E-and F (also from left to right). Plots C and 
D are being adversely affected by vignetting at the edge of the 70 mm transparency. Original negative scale - 1:6,330; 
north is at top of page. 
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“ Feasibility of Detecting Major Air Pollutants 
by Earth-Oriented Satellite-Borne Sensors ,* 1 


Investigators: Dr. Harold J. Paulus , David W. Hoffman, 

Department of Environmental Health 


To further aid in the evaluation of air pollution in the Minne- 
sota area, ERTS-1 MSS imagery are being sought for areas where enor- 
mous sources of air pollution are present. The identification and 
densitometric examination of air pollutants, especially suspended 
particulates, should help to evaluate the lesser sources of air 
pollution in Minnesota. 

Contrast reductions over the metropolitan areas are being 
examined for possible construction of isopleths of particulate 
matter, and for correlation with air pollution index values fron 
local measurements. 

Data sources fron which atmospheric parameters are directly 
available are the Minnesota Pollution Control Agency, Minneapolis 

\ 

Air Pollution Control Division, St. Paul Air Pollution Control Division 
and local meteorological stations. 

The examination of pollution injury to vegetation (crops and 
forests) in the area of air pollution sources will be evaluated. 

The capabilities of the ERTS-1 MSS imagery relative to air 
pollution will soon be presented to the Minnesota Pollution Control 
Agency and other local agencies concerned with air pollution. 
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USE OP REMOTE SENSING IN THE ANALYSIS OP SNOWMELT FLOODS 

by 

C. Edward Bowers 

I Introduction 

The objective of this study is the evaluation of remote sensing data 
and particularly ERTS Imagery as an aid to the determination of percent snow 
cover and possibly water content of snow on the ground. This information Is 
of interest in the application of some mathematical simulation models to the 
forecasting of spring snowmelt floods. 

The largest floods of record for most areas in Minnesota have usually 
occurred in the spring and consist in whole or in part of snow melt runoff. 

r 

Thus, in attempts to predict the magnitude and timing of the spring floods 
it is necessary to obtain data on the snowpack during late winter and early 
spring. If small areas are of interest, it is relatively easy to make 'field 
measurements of water content and estimates of percent snow cover. However, 
the primary problems usually develop in the 16,200 square mile Minnesota 
River Basin, the 19,100 square mile area of the Mississippi River above 
Anoka or combinations of these two areas and a downstream segment for a total 
of 59,200 square miles. A similar problem exists for 40,000 square miles of 
the Red River of the North Basin. Thus, data collection for areas of this mag- 
nitude can be formidable. 

Currently, field sampling of the water content of snow at 60 locations 
in Minnesota, North Dakota and Wisconsin provide some data on water content 
of snow. (This corresponds to 1 measurement for about 1500 square miles.) 
Information on percent snow cover is not generally available. 
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II - Current Research 

Studies of snowmelt flood analysis in the upper midwest are currently 

underway at the St. Anthony Falls Hydraulic Laboratory, Department of Civil 
and Mineral Engineering, University of Minnesota. The studies sponsored by 
OWRR involve fitting of a mathematical simulation model to the Minnesota 
River Basin, using available hydrologic data for past years. The model, 
involving about 8,000 input statements, requires information on water content 
of snow and percent snow cover as well as river discharge, precipitation, 
temperature and related variables. Snow melt can be computed on a continuous 
basis either by a degree-day or temperature method, or an energy basis involv- 
ing solar radiation as well as temperature. The model computes the melt and 
the balance of snow in storage but must be provided with a relationship between 
water content of snow and percent snow cover. Figure 2 illustrates one simple 
relationship; this assumes that if the snow water-content is greater than 1 
inch, there is 100$ snow cover, with a linear relationship between 0 and 1 . 
inch. This has not been verified and is the area in which remote sensing, 
applied to the 16,200 square mile basin should be helpful. Figure 3> shows 

part of the print-out of the model and Figures 4 and 5 show computed or 
* 

simulated discharge compared to observed discharge. The results for the 
Pomrne de Terre River (1967) Figure 4, are good. The agreement for the 
Minnesota River at Jordan (Figure 4) and the Chippewa River (1969), Figure 5* 
indicate the need for improvement. 

It is hoped that the ERTS Imagery for the Minnesota Basin in late winter 
and early spring will provide a sound basis for estimating percent snow cover 
over the large areas involved in the mathematical simulation study. . However, 
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there is no reason to limit the ERTS study to the Minnesota Basin. Since 
field data on water content are available for all of Minnesota and parts of 
Wisconsin arid North Dakota, these additional areas will also be considered 
and possibly used in the study. 

Flood damages in Minnesota amounted to $92,000,000 in 19^5 and $77>000,000 
in 1969, two of our worst flood years. Annual flood damages are on the order 
of 10 to 15 million dollars. Research studies plus flood plain management 
should help to reduce these figures. 
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Typical Output from SSARR model, 
Chippewa River Watershed 
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Fig, 4. Simulated and Observed 
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Flow at Two Stations in Minnesota Paver Basin 
in 1967. 
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11 A Proposal for Precambrian Mineral Resource 
Evaluation in Minnesota utilizing ERTS Imagery . 

Investigator: Dr. Donald Davidson, Jr. 

Department of Geology 

I propose an investigation of the discriminatory ability of ERTS 
imagery to geological relationships in Precambrian rocks of Minnesota. 

These rock units hold promise for potential mineral resources as do 
important ore deposits found in similar rocks to the north in Canada. 

Two major Precambrian rock units are proposed for examination using 
ERTS imagery. The research planned involves the discrimination of rock 
types to show their aerial extent and an interpretation of the structural 
relationships between and within the various rock units. A considerable 
amount of ground truth information developed by the bedrock mapping program 
done in the past 10 years by the Minnesota Geological Survey is available 
for both units proposed. 

The two major rock units to be examined are: 

1. Archean Greenstone Belts 

Greenstone belts throughout Canada are the most important sources of 
Cu-Ni-Zn-Ag-An ores in that country. Several greenstone belts are known to 
occur in Northern Minnesota. The best studied belt, the Vermilion belt, has 
proven to be essentially barren. Moreover, geologic mapping and exploration 
work has been conducted in most known greenstone belts. However, adequate 
study and evaluation of these belts has been hampered by the lack of regional 
as well as detailed information concerning distribution, structural relation- 
ships and other geological data. Regional structural and lithologic analysis 
In and near the several belts may hold the key to either why a belt is 
barren or where more favorable areas of mineralization are to be found. 
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2. The Troctolite Series of the Duluth Gabbro Complex 
Known Cu-Ni mineralization is found in the lower, western part of 
the Duluth Gabbro Complex in rocks of troctolitic composition. To date 
Cu-Ni ore bodies are found in the basal contact zones of within the 
Troctolite Series. The total extent of the boundaries of the troctolitic 
units have not been determined primarily because of their extensive 
occurrence. Since the basal contact of the units is where known ore deposits 
have been found, regional patterns, structural studies and even a broad 
delineation of the major rock units could be of great value to continued 
exploration. Additional information of regional scope could greatly enhance 
the chance for further success in discovery of Cu-Ni deposits. 

This proposed program, while having considerable scientific value, 
merits consideration because the scientific results can be expected to 
have considerable value in determining the metal resource potential of the 
Duluth Gabbro Complex. \ 



